Parity oscillations of Kondo temperature in a single molecule break junction 
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We study the Kondo temperature (Tk ) of a single molecule break junction. By employing a numerical renor- 
malization group calculations we have found that depends dramatically upon the position of the molecule in 
the wire formed between the contacts. We show that Tk exhibits strong oscillations when the parity of the left 
and/or right number of atomic sites (N^, Ng) is changed. For a given set of parameters, the maximum value of 
Tk occurs for (odd, odd) combination, while its minimum values is observed for (even, even). These oscillations 
are fully understood in terms of the effective hybridization function. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.21.Ac, 73.21.Hb, 73.21.La, 73.63.Kv, 73.63.Nm, 73.21.La 



Kondo effect (KE) is one of the most intriguing phenom- 
ena of strong correlated systems,' which was beautifully ex- 
plained by J. Kondo in the 60's in the seminal theoretical 
work on the minimal resistance in magnetic alloys.— KE has 
revived in the later 90's with the advent of the scanning tun- 
neling microscope (STM) that has facilitated the manipulation 
of the matter at atomic scale. For instance, STM has allowed 
observation of interesting facets of the KE in quantum dots 
(QD)^ and in single atom or molecule on metallic surfaces^— 
which have motivated a huge number of experimental^ and 
theoreticalSi^ investigations. 

One of the experimentally accessible signatures of the KE 
in nanoscopic system such as QD or magnetic molecules at- 
tached to metallic contacts is the strong modification in the 
conductance across the system, observable when the system 
is cooled down below the so-called Kondo temperature (Tk)- 
In QD, for instance, Tk is found to be in the sub Kelvin re- 
gion whereas for large molecules attached to metallic leads 
Tk can be much larger)^ In both cases, in the limit of very 
strong Coulomb interaction, Tk depends strongly upon the ef- 
fective hybridization (A) that connects the localized magnetic 
moments to the conduction electrons^JS asA Tk ~ exp(7TEd/ A), 
where Sd (< 0) is the energy of the localized orbital respect to 
the Fermi level. Controlling A or Sd is therefore crucial for ob- 
taining higher Tk, which is fundamental for possible techno- 
logical application of KE. While tuning Sd is relatively simple 
in QDs by mean of gate voltages, in molecules, on the other 
hand, it becomes a more complicated task. Conversely, geo- 
metrical parameters are more suitably modified in molecules 
than in QDs and has proven to produce important modifica- 
tions in Tk via hybridization functionJi 

A suitable experimental arrangement to study the KE is the 
break junction (BJ) molecular structures, in which a metallic 
wire (gold wire, for instance) is stretched until a few-atom ID 
chain bridges the gap between the electrodes before the com- 
plete break up of the wire; Owing to the dependence of Tk 
upon A it has been shown that Tk can be mechanically mod- 
ulated in BJ experiments^' by changing the distance between 
the electrodes. Motivated by this experiment, in the present 
work we study the Kondo temperature of a spin- 1/2 magnetic 
impurity coupled to metallic contacts through two finite (left 
and right) quantum wires (QW), as illustrated in Fig. ^a). By 
employing a numerical renormalization grou p '^i'^ (NRG) cal- 
culation we find a strong dependence of Tk upon the parity 



of the number of sites (Nr,Nl) of each QW as well as their 
length. The dependence upon the (Nl,Nr) parity combina- 
tion results in an oscillating behavior of Tk as function of Nl 
or Nr, akin to what has been observed in Manganese phthalo- 
cyanine (MnPc) molecules on top of Pb islands, reported in 
Ref.[l^ Although the system under investigation here is quite 
different from the one studied in Ref. 19, the origin of the os- 
cillation of Tk can be interpreted likewise. While in their case 
the enhancement of A originates from the formation of mul- 
tiple quantum well states between the Pb atomic layers, here 
the enhancement of A results from the localized states of the 
atomic sites of the QW. 
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FIG. 1. (Color online) (a) Illustration of a Au quantum wire coupled 
to metallic contacts with a embedded Cgo molecule, (b) Pictorial 
representation of the model. The site labeled as "0" represents the 
impurity site with strong on-site coulomb repulsion. 



Our system model is schematically represented in Fig. [Hb) 
and is modeled by the Anderson-type Hamiltonian that can be 
split into five terms as 



(1) 



where Himp, Hcb and Hmres describe, respectively, the inter- 
acting impurity, the free electrons in the conduction bands and 



the electrons in the wires, //„„y,-u/>f.s couples the impurity to 
the two wires and Hch-mres couples the wires the their respec- 
tive conduction bands. In terms of creation and annihilation 
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fermion operators the Hamiltonians read 
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In Eqs. |2]|6l the operators c^^ (Qo-) creates (annihilates) an 
electron in the orbital d with energy Sd, c].^^ {ctka) creates 
(annihilates) an electron in the fth conduction band with en- 
ergy et, c'.^^ (cifo-) creates (annihilates) and electron in the ith 
site of the fth QW with energy eo spin cr. Finally, t is the 
hopping between two adjacent sites in the wires and Vrk and t' 
couple the QWs to their conduction bands and to the impurity, 
respectively. The conduction bands are characterized by a flat 
density of states, picij) = (1/2D)0(D - where D is their 
half bandwidth and 0(x) is the Heaviside step function. It is 
worth emphasizing that although the motivating experiment 
was realized using Ceo molecule coupled to Au metallic con- 
tacts, this model is rather general. In the particular context of 
molecular BJ, vibrations may be important in certain range of 
parameters, but this aspect is beyond the scope of the present 
work. 

In order to properly address the Kondo physics of the sys- 
tem, the full Hamiltonian is approached by using the nu- 
merical renormaUzation method, with which we can calcu- 
late the thermodynamical properties. Within the NRG ap- 
proach we discretize the effective conduction band "seen" 
by the interacting impurity. The effective conduction band 
can be determined by exact calculation of the local non- 
interacting {U - 0) Green's function (suppressing the spin 
index), gdd{(^) - - + ^'^)] where 2 = 2fi((y) + SL((y), 
in which the fth self-energy is given by 



(7) 
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being the diagonal GF associated to the unperturbed conduc- 
tion electrons in the leads. The fraction in Eq.|7]is continued 
until all the sites of the of the fth wire and the ^th conduction 
band are taken into account. 



The hybridization of the localized orbital 'W with the ef- 
fective band is give by A(a)) = Im[2(w)] = A^(cli) + A«(w) 
[Hereafter we will refers to A(0) just as A]. The hybridization 
function is logarithmically discretize d'^'^" to map the system 
in a Wilson's chain form^^ 

DO 

H = Hi,„p + t'Yj i'^do-'^Oa- + C(,^Crfo-) + ^ Eio-cj^Cia- 
cr i=o 

+Xi^'(^''-'''+''"+''''+i<^'''^)' 



where f,'s are calculated via A(w), following the recipes 
described in Ref. Once we have mapped the system 
in the Wilson's form, we proceed the NRG calculation, 
which is based in the iterative diagonalization of the effective 
Hamiltonian. "* After reaching the strong coupling fixed point 
we can calculate the magnetic moment within the canonical 
ensemble as 



^ ^ [(v|52|v) - (v\S:,\y)f] e-^"""\ (10) 



where is the Boltzmann constant, Z(T) - 

YjvSxp{-Ey/kBT) is the canonical partition function, 
is the spin operator, |v) and Ey are, respectively, the 
eigenvector and its corresponding eigenvalue of the full 
interacting Hamiltonian, which are naturally calculated 
in the NRG procedure.— Following Wilson's criterion, 
we define Tk from the "impurity" magnetic moment as 
kBTKXimpiTK) - 0. 07 07 (gfiB)^, (that is the magnetic moment 
of the full system subtracted by the contribution of the 
effective conduction band), g is the electron ^-factor and fig 
is the Bohr magneton. 

Before starting the presentation of our numerical results, 
lets analyze the hybridization function at the Fermi level, 
which is the most relevant parameter to determine the behav- 
ior of T/f in our calculations. It is straightforward to show 
from Eq. Qthat A possesses only three distinct vales, 

Amin - Aq for (even,even) 

ii„i = Ao (5- + ff) for (even,odd) or (odd,even) , 
A„,flj = 2aAo for(odd,odd) 

(11) 

where we have defined Ao - nt'^/D and denoted A„„„, A„„ and 
Amax, the minimum, intermediate and maximum value of A, 
respectively, and a = 2{D/nt)^ is a dimensionless parameter 
that can be modified, for instance, by stretching the QW as is 
was done in the Ref. [U To obtain our numerical results lets 
set D = 1 as our energy scale. With that we choose hereafter 
(unless otherwise stated) U = 0.5, £d - -0.25, eo = (at 
the particle-hole symmetric point), Vr - Vl - t - 0.15 and 
f = 0.1. 

In Fig.|2]we show the hybridization function vs energy for 
various values of Nl and Nr. In Figs.|2la),|2|b), and|2{c) we 
fix = and show A(w) for Nr - I, Nr - 2 and Nr - 3, 
while in Figs.|2td),|2|e), and|2tf) we keep Nl - I fixed and 
show A{oj) for Nr - 2, Nr - 3 and Nr - 4. The number of 
peaks of A{(d) is given by maxiNi, Nr) for equal parity and 
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FIG. 2. (Color online) Hybridization function vs energy for various 
values of Ni and A'^ [denoted in the figure as (Ni, Nr) using = 
= f = 0. 15 and t' = 0.1. Notice that A possesses three difi"erent 
values, depending on the parity of Nl and Nr. The minimum (A,,,,,,) 
and maximum (A,„a j) value of A is obtained for (even, even) (b) and 
(odd, odd) (e), respectively, while for all the other combinations A 
has an intermediate value, A„„. 
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FIG. 3. (Color online) Magnetic moment as function of temperature 
for various values of = Ni = N. Panel (a) and (b) correspond to 
N even and odd, respectively (see values of A' in the legends). 



Nr + Nl for different parities. Although the structure of A(w) 
away from the Fermi level has some effect on Tk, the most rel- 
evant contribution comes from the structures at or very close 
the the Fermi level. For the parameters set above, we obtain 
A„„„ ^ 0.0314, A„,„., ^ 0.566 and A,„, ^ 0.299. These distinct 
values of A are crucial for determining the Kondo temperature 
of the system, which in our case can be roughly estimated as-- 
ksTKiA) ~ exp[-7rL'/8A]. It is clear that Tk increases as A in- 
creases. 

In Fig. [3ja) and Fig. [3jb) we show the magnetic moment 
as function temperature for various values of N = Nl = Nr 
(the symmetric case) even and odd, respectively. The case of 
N - [Fig. Oa), O (black)] corresponds to the single im- 
purity coupled to two conduction band. The low temperature 
suppression in the magnetic moment results from the Kondo 
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FIG. 4. (Color online) Kondo temperature as function of number of 
sites (Nr) for a fixed Ni = I. O (black) and □ (red) coiTespond 
to even (odd) Nr, respectively. The zig-zag (green) line shows the 
(even, odd) oscillation, similar to those observed in Ref.lTgi 



screening of the local spin [these curves are used to Tk, as 
discussed above]. On the other hand, in the high temperature 
limit the ksTx — > (gfiB)'/^, as expected. Notice in Fig.[3ta) 
that Tfc increases when (even) increases. Conversely, Tk 
decreases when (odd) increases as seen in Fig. Ob). No- 
tice also that Tk can be at least two order of magnitude larger 
for A^ odd than for A^ even. This huge difference will be an- 
alyze below. The result for A^ = 1 [Fig. [3jb), O (black)] is 
equivalent to those reported in Ref.[l3 The negative values of 
ksTximp within a small range of T results from the subtraction 
of the effective conduction band contribution. 

In order to show the behavior of the Kondo temperature for 
larger and different values of Nl and Nr we show in Fig. |4] 
Tk as function of Nr for a fixed number Nl even (left) and 
Nl odd (right). The O (black) curves correspond to A^r even, 
while □ (red) curves corresponds to Nr odd. Notice that for 
Nl even [(Fig.|4ja),|4|b) and|4jc)] Tk increases with Nr even 
[O (black)] while it decreases for Nr odd [□ (red)]. Observe 
again that for small Nr Tk is almost two order of magnitude 
larger for Nr odd than for Nr even (keeping Nl even). This 
difference decreases asymptotically for large Nr and vanishes 
asymptotically as A^ oo. This results from the fact that 
in this situation the conduction electrons near the Fermi level 
are more strongly coupled to the impurity, reflecting the fact 
that A,„, is larger than A„„„ as clearly shown in Fig. |2] For 
A^L odd (Fig.|2jd),|2te) and|2tf) we observe a similar behavior 
(Tk increases as Nr even increases and decreases as Nr even 
increases) but in this case the curves collapse onto each other 
very quickly (typically for A^^ = 10) to a larger value, when 
compared to the case of A^^ even. At least for small Nl and Nr 
we can roughly estimate the ratio between Tk^ for the three 
distinct values of A as 



TKiAb) 



(12) 



where a and b stand for min, int and max. Using the param- 
eters chosen above we obtain TK(Amt)ITK(Amm) ~ 2.7 x 10", 
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FIG. 5. (Color online) Kondo temperature as function of £,/ (eq = 0) 
(a) and eo (srf = -U 12) (b) for various configuration of {Nl,Nr) as 
shown in the legend. The other parameter are the same as in the 
previous figures. 



while TK(A,„ax)/TK{Ai„,) ^ 1.36. These values are consistent 
with the huge difference between the values shown in □ (red) 
and O (black) curves of Figs.|4la),|4lb) andHJc) and small dif- 
ference in the related curves of Figs. EId),|4le) and^f). The 
behavior of Tk with increasing Nl and Nr for the same par- 
ity combination cannot be explained in terms of A. This can 
be reasonably understood in terms of the formation of a small 
sub-band inside the conduction band, due to a large number 
of atomic sites and the energy dependence of hybridization 
function near the Fermi level. In the Umit of Nl,Nr oo the 
sub-band becomes a smooth curve, resulting in a indepen- 
dent of the lengths of the wires. The zig-zag (green) line in 
Fig. |4jb) shows the even-odd oscillations in Tj^, very similar 
to what was observed in Ref. 19. Finally, in Fig. |5]we show 
the robustness of these results against particle-hole symmetry 
breaking. In Fig. |3a) show Tk as function of Sd for eo - 0. 
Notice that, although more pronounced for the (even, even) 
case, Tk increases as is shifted upward from -U/2 for all 



parities [(0,2), (1,2) and (1,3)]. Same behavior is obtained 
for the other side (not shown). These are consistent with the 
general expression,^-' Tk ~ Exp[-7r\sd\{sd + U)/(2AU)] for 
constant hybridization function. When we keep Sd and vary 
£() about the Fermi level [Fig. ISjb)] we see that Tk increases 
for (0,2) but decrease sUghtly for (1,2) and (1,3). This re- 
sults from the fact that, as eo deviates from the Fermi level, 
A(0) decreases if A{a>) possesses a peak at the Fermi level 
as in the (even, odd), (odd, even) or (odd, odd) cases, but in- 
creases when A{a>) exhibits a dip at the Fermi level as in the 
(even, even) configuration. 

In conclusion, we have presented a detailed study of the 
Kondo temperature of a single molecule break junction. By 
employing a numerical renormalization group we show that 
Tk is strongly dependent of the parity of the number of atomic 
sites in each piece of QW connecting the molecule to the con- 
tacts. More interesting, we show that the Tk oscillates when 
the parity of the number of sites of the wires changes. These 
oscillations are interpreted in terms of the effective hybridiza- 
tion function A(w). For (even, even) and (odd, odd) configu- 
rations the effective coupling A is minimum and maximum, 
respectively, while for (even, odd) or (odd, even) configura- 
tions A possesses an intermediate value. Within this picture, 
the huge variation of Tk is readily estimated by a simple an- 
alytical calculation, which can vary up to a factor of 10^ [in 
the case of changing from (even, even) to (even, odd)]. Our 
results provide a very clear picture of the main ingredient re- 
sponsible for the dramatic dependence of Tk on geometrical 
configuration of single molecule break junctions as well as of 
magnetic molecule on atomic layer surfaces. Moreover, we 
believe our results can be used to guide experimental realiza- 
tions of high-T^: experiments. 
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